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Abstract 

Second initial boundary problem in narrow domains of width e ^ 1 for 
linear second order differential equations with nonlinear boundary condi- 
^ ' tions is considered in this paper. Using probabilistic methods we show 

•^O , that the solution of such a problem converges as e 4- to the solution of 

a standard reaction-diffusion equation in a domain of reduced dimension. 
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rvj ■ This reduction allows to obtain some results concerning wave front prop- 
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agation in narrow domains. In particular, we describe conditions leading 
("^ ' to jumps of the wave front. 

>--^ ■ Keywords: Reaction Diffusion equations, Narrow Domains, Wave Front 

Propagation, Instantaneous Reflection 
> . 



This is an electronic reprint of the original article published by the 



Asymptotic Analysis Volume 59, Number 3-4 / 2008, pp.227-249| This 



reprint differs from the original in pagination, typographic detail and 
typo corrections. 



2 


in (0,T) X D^ 


/(^), 


on {0} X D' 


-ec(x,y,'u')u', 


on (0, T) X dD' 



M.Freidlin and K. Spiliopoulos / RDE's in Narrow Domains and Wave Fronts 

1 Introduction 

For each x £ M", let D^ be a bounded domain in R™ with a smooth boundary 
dDx. Assume, for brevity, that Dx is homeomorphic to a baU in M™ and 
contains € M™. Consider the domain D = {{x,y) : a; € R",?/ G ZJ^} C K"+™. 
Assume that the boundary dD of D is smooth enough and denote by 7(0;, y) the 
inward unit normal to dD. Assume that "f(x,y) is not parallel to the subspace 
R" c R"+'" for any {x,y)edD. 

Denote by 13'^, < e << 1, the domain in R"+'" obtained from D by con- 
traction: D'^ = {{x, y) : X G R", ye~^ G Dj.}. If n = 1, D'^ is a narrow tube (or 
a strip for m = 1) for < e << 1. If n > 1, then D'^ is a thin layer. 

Consider the problem: 

ut = ^Au', in(0,T)xD^ (1) 

u'{0,x,y) 

where 7*^ is the inward unit normal to dD'^ . The functions / and c are sufficiently 
regular and bounded; / is assumed to be nonnegative. Our goal in this paper 
is to study the behavior of solution of problem ([1]) as e | 0. Using probabilistic 
methods, we will prove that u'^(t,x,y) converges as e J, to the solution u(t,x) 
of the problem: 

Ut = l-A,u+l-V{logV{x))V,u+l§^c{x,0,u)u, in(0,T)xR" 
2 2 2 V (x) 

u{0,x) = fix), on{0}xR". (2) 

Here V{x) is the volume of D^ and S{x) is the surface area of dD^. One 
can expect that, under certain assumptions on the nonlinear term c(x,y,u)u 
in ([l}, the solution u'^(t,x,y) can be approximated by a running- wave- type 
solution. Corresponding results on the standard reaction diffusion equation ([2]) 
(see chapter 6 and 7 in 0]) allow to describe the asymptotic wavefront motion 
for ([T]). We will see how the motion of the interface (wavefront) depends on 
the behavior of the cross-sections D^ of the domain D. In particular, using the 
results of 4, (chapter 6) we will see that in the case of the nonlinear term of 
K-P-P type the wavefront can have jumps. 

Consider the Wiener process {X^,Y^'^) in D'^ with instantaneous normal re- 
flection on the boundary of D'^. Its trajectories can be described by the stochas- 
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tic differential equations: 

X^ = x + W}+ j Yi{Xl,Y:)dLl 
Jo 

y/ = y + W?+ f -,',{Xl,Y:)dLl. (3) 

Jo 

Here W^ and W^ are independent Wiener process in R" and M™ respectively 
and (a;, y) is a point inside D'^. Moreover 7I and 7I are projections of the unit 
inward normal vector to 9C on R" and R™ respectively. It is easy to see that 
lime^o |e~^7f| = -prr and linie^o I72I — 1; where | • | denotes Euclidean length. 
Furthermore LI is the local time for the process (Xj,F/) on dD'^, i.e. it is 
a continuous, non-decreasing process that increases only when {X^^Y^'^) G dD^ 
such that the Lebesque measure A{t > : {X^, Y^) e dD'^} = (see for instance 

If {Xf,Y^'^) is defined by ^, then as it can be derived from Theorem 2.5.1 
in [3], u'^{t, X, y) satisfies the following integral equation in the functional space: 

u'{t,x,y)^E.,,yf{Xt)exp[f ec(X|, F/, zi'(i - s,X|, r/))dLJ], (4) 

Jo 

where E^^y denotes expectation and the subscript {x, y) denotes the initial point 
of {XI, F/). Equation ^ has a unique solution if, say, c{x, y, u) has a bounded 
derivative in u. 

Let Xt be the solution of the stochastic differential equation 

Xt^x + W} + j U/ {\ogV {X,))da. (5) 

Then the solution u{t,x) of equation ^ satisfies the equality: 

u(t,x)^E^,f{Xt)exp[J^ lp^c{X,,0,u{t-s,Xs))ds]. (6) 

We prove that the component X^ of the process {X^,Y^) converges in a 
certain sense to Xf. This together with uniform in < e < 1 bounds for 
u'(t, X, y) and its derivatives allow to prove that the solution of ^ converges to 
the solution of ([61) as e | uniformly on each compact subset of [0, 00) x R"+"'. 

In the next section we consider averaging of integrals in local time. This 
result allows in section 3 to prove convergence of the integral in the right side 
of the first of equations in ^ to the integral term in ([5|) and convergence of 
exponents in Q and ^. Together with a-priori bounds obtained in section 
3, this implies convergence of u'^{t,x,y) to u{t,x). Some results concerning 
wavefront propagation are presented in section 4. 
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2 Averaging of Integrals in Local Time 

Let H{x^ y) be a smooth and bounded function. We want to consider the hmiting 
behavior as e | of expressions hke /„ eH{Xl, Y^ /e)dL\ (see Lemma 2.1 below). 
We will assume that the unit inward normal ^{x,y) to dD and the function 
H{x, y) are both three times differentiable in x and y. 

Lemma 2.1. Define Q{x) = yj^ Jgj^ H{x, y)dSx, where dSx is the surface 
element on dD^- Then for every T > and small enough e, there exists a 
constant K independent of e such that: 

(i). snpo<t<T E\ /o ^Q{XI)ds - /J ei/(X|, r//6)|7|(X|, r/)|dL||2 < Ke\ 
(ii). For every ^ > we have 

^{supo<,<T I /o hQi^Dds^Io 6H(X|,y//e)|7|(X|,y/)|dL|| ><5} < K^. 
The proof of lemma 2.1 relies on the following lemma, which we prove first. 

Lemma 2.2. For every T > and small enough e, there exists a constant 
Ki independent of e such that: 

EleLfpl"^ < Ki 

Proof. 

Consider the auxiliary problem 

^yv{x,y) = Q(a-), y e I?, C M™ 
dyv{x,y) 



dn{x,y) 

v) 
\ili^,y)\ 



1, y^dDx, (7) 



where nix, y) — ?tf'^L and x € M" is a parameter. Let 



where S{x) is the surface area of D^ and V{x) is the volume of D^. As it can 
be derived from [T], a smooth in x and y solution v{x,y) of problem ([7]) exists 
and is bounded together with its first and second derivatives. So we can apply 
Ito formula to the function ev{x,y/e), and get: 

eMX^,Y,ye) = eMx,y/e) + j ^\^MXl,Y:ie)ds + j \^yv{Xl,Y: /e)ds 

+ / ^{Vxv{XlY:it), dWl) + / e{Vyv{Xl Y^/e), dW^) 
Jo Jo 
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e\V.v{Xl, y//e), 7i^(X|, Y:))dLl 



e{yyv{Xl,Y^/e), 7^(X|, Y,^))dL\ (9) 

Recalling now that linij^o k^^7il = ttt a-nd linie^o I72I — 1 ^.nd that v satisfies 
([7]) one easily concludes that there is an eg — ^q{\\\^ xv\\\t^\) > such that for 
every e < eo- 

+ e^\Vyv\'\\T+\\^QfT], 

where for any function g, \\g\\ = sup^ l5(-z)|- Here, we also used the fact that the 
local time is increasing function of t. This proves the statement of the lemma. 

D 

Proof of Lemma 2.1. We consider the auxiliary problem 



Ayv{x,y) = Q{x), y e D,r C 
dyv{x,y) _ 



dn{x,y) 



H{x,y), yedD,, (10) 



where n(x, v) — ,^1, '^^1 and x € M" is a parameter. 

The necessary and sufficient condition for the existence of a solution for (jlOp 
is that 

Q^'^^^TTT-^f H{x,y)dS^, (11) 

where dSx is the surface element on dDx and V{x) — volume(I?a:)- 

Applying Ito formula to the function ev(x, y/e) and using the bounds ob- 
tained in Lemma 2.2 wc get the following inequalities: 

sup E\ f ^Q{Xt)ds- f eH{Xt,Y:/eMiXt,Y:)\dLlf< 

0<t<T Jo ^ Jo 

< e^C{2\\v^ + W^AxvfT^ + \\\Vxv\'\\T + \\\Vxv\\\'Ki) + e^C\\\Vyv\^\\T, 

which proves statement (i) of the lemma. 

For part (ii) one makes use of the Doob maximal inequalities (see [TT] , page 
14): 



E[ sup I / iVMXt,Y:/e),dW, 

0<t<T Jo 

E[ sup I [ {Wyv{XlY:/e),dW, 

Q<t<T Jo 



2^112 < 4 II |Vj,wp II T 
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Then, following the procedure that proved part (i) we get that there is an eg > 
such that for every e < eo: 

E[ sup I / ^QiXDds^ f eH{XI,Y,yeMiXt,Y,^)\dLl\r < 

0<t<T Jo ^ Jo 

< e''Ci4\\v^ + \\^A,v\\'T'+A\\\V,v\^\\T+\\\V,v\fK,) + e^C\\\\/yv\^\\T, 
which together with Chebyshev inequality proves statement (ii) of the lemma. 

n 



3 Limit of u'^. 

In this section we consider the limit as e — >■ of the solution u*^ to problem ^ . 
The result is given in Theorem 3.4. The proof will proceed as follows. First (in 
Proposition 3.2.) we write down an integral equation in the space of trajectories 
for the solution of ([T]). Then in Lemma 3.3 we consider the mean square limit as 
e — > of the underlying stochastic process with instantaneous normal reflection 
on the boundary of D"^ (see ([3])). Lastly an important ingredient to the proof 
are the a-priori bounds for u' and its derivatives. These a-priori bounds are 
independent of e, their derivation is standard and are given for completeness in 
Proposition 3.7. 

We assume that the initial function f{x) of problem ([1]) is bounded, non- 
negative and can have finite number of simple discontinuities. The function 
c(a;, y, u) is assumed to be uniformly bounded in all arguments, continuous in 
x,y, Lipschitz continuous in u and that there exist constants M, A^ > such 
that c{-,-,u) < -M for u> N. 

In addition we assume that the boundary of D^ satisfies dD^ E C''+°(M'"), 
where a £ (0, 1). 

Remark 3.1. For the existence of a classical solution to ([T]) one actually 
needs only to assume dD^ e e2+"(]R"). The assumption dD^ e e^+''{R"') is 
being done solely for the purpose of Lemma 3.2 and Theorem 3.3. 

Let (X%y%L") in R" x W" x R\_ satisfy ©. Then we have: 

Proposition 3.2. Problem ([1]) has a unique classical solution in [0,T) x D'^ 
which satisfies: 



u'{t, X, y) = E,^yf{Xt) exp[ / ec(X|, Y^,u\t - s, XI, Y!))dL\ 

Jo 



(12) 
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Proof. Under our assumptions, the uniqueness and existence of a classical 
solution to (m follows from Theorem 7.5.13 of [^). The equation (fT2)) follows 
from Theorem 2.5.1 of j4j. 

n 

In order now to consider the limit as e — )> of (|12p . we need first to examine 
the asymptotic behavior of X^ as e — t- 0. 

We will prove that X^ converges as e 4- to Xj , where Xt is the solution to 

/■* 1 
Xt^x + Wl+ -V(log V{X,))ds, (13) 

Jo ^ 

where V{x) — volume(Dj;). Hence, we see that as e J, 0, the effect of the 
reflection on the boundary is an extra drift term. A sketch of the proof for the 
above result is given in chapter 7 of 6 . More details are given here. 

Lemma 3.3. For any T > we have 

sup E^\Xl - Xt'i^ ^ Q as e^O. (14) 

0<t<T 

Proof. It is not difficult to see that 75^(x,2/) = g /ifa'^^i ItIC^i V)\ ^^id that 



9Dx 172 



^i^'^'^^d^, = Vy(x). (15) 



\l\{x^V)\ 



Then Lemma 2.1 with H{x,y) = i'^i('^\i and Q{x) — Vlogl/(x) implies that 



for small enough e there exists a constant K independent of e such that 
sup E\ f lviog{ViXt))ds~ f jl{Xl,Y^^)dLl\'<e'K. 

<t<T Jo ^ Jo 



0<t<T 

Now we write 



(16) 



X^ -Xt ^ J lliXl, Y:)dLl - / ^ V \og{V{X,))ds 
^ [J^ jliXl, Y:)dLl - J^ iv log{V{Xl))ds] 



' iv log{V{x:))ds - / Jv log{V{X,))ds] (17) 
^ Jo ^ 



Then Gronwall Lemma and (|16l) give: 



sup E^\X^ - XtY ^0 as e ^ 0, (18) 

0<t<T 



which is the statement of the lemma. 
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n 



Consider now the solution, u, to the equation 



u{t,x)=E,fiXt)e^p[[ c{X,,u{t-s,X,))ds], (19) 

Jo 

where 

cix,u{t,x)) = lp^c{x,0,u{t,x)). (20) 

For notational convenience we will also denote c{t, x) — c{x, u(t, x)). 

Since c{x, u) is Lipschitz continuous in u, the solution of (1191) exists and is 
unique. Our assumptions on the functions /, c and the boundary dDx, imply 
that the solution u to ([T9| is actually the classical solution of the following 
parabolic problem: 

ut ^ l-A^u + l-V{logV{x)) ■ V^u + l-§^c{x,0,u{t,x))u, in(0,r)xM" 
2 2 2 V (x) 

u{0,x)^f{x), on{0}xM". (21) 



Theorem 3.4. Under our assumptions, we have that 

u'^{t, X, y) —> u{t, x) as e — ?► 0, uniformly in any compact sunset of M-|_ xM" xR™, 

where u'^ {t , x , y) , u{t,x) are the solutions to dD) and (ICT respectively. 

Proof. By Proposition 3.6 and the well known theorem of Ascoli-Arzela 
we get that there exists a subsequence of {u'} (which for convenience we will 
denote again by {u*^}) and a function u, such that: 

u*^ — > u as e — >■ 0, uniformly in compacts. 

We will prove that u actually satisfies ([T^ which then implies that u satisfies 
(|2T|) . Fix t and x and consider the solution v{y) — v'^'*-'^{y) to the elliptic 
boundary value problem: 

Ayv{y) = c%t,x), yeD^cR"" 
dyv{y) _ 1 



dn{x,y) \ll{x,<iy)\ 

Problem (|22]) is solvable if 



c{x, ey, u'-{t, X, ey)), y G dD^. (22) 



1 /■ 1 

c'(i,x) = — r^ / — -c{x,ey,u'{t,x,ey))dS.:,. (23) 
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Proceeding similarly now to Lemma 2.1 and recalling that v satisfies (|22p . we 
see that there is a constant X' = K{\\V^v'\\,\\VyV^\\,\\A^v''\\,\\Ayv'\\,\\vl\\,\\-fl\\,T) 
such that: 

sup E\ [ y{t-s,Xt)ds- [ ec{X:,Y,\u'{t - s,Xt,Y:))dLlf 

Q<t<T Jo ^ Jo 

< e^K'{l+ sup E[eLlf) (24) 

0<t<T 

We observe that K"^ depends on e only through functions that are uniformly 
bounded in e (Proposition 3.7). This observation and Lemma 2.2 imply that as 
e-)-0: 

sup E\ [ y{t-s,Xt)ds- [ ec{Xt,Y:,u'{t-s,Xt,Y,'))dLlf^O. (25) 

0<t<T Jo ^ Jo 

Moreover the Lebesque dominated convergence Theorem, Lemma 3.3, the com- 
pactness of the family {u^} and (1^ . imply that as e — > 0: 

sup E\ [ ^c^{t-s,XI)ds- I c{t-s,Xs)ds\^ ^0, (26) 

0<t<T Jo 2 Jo 

where c'^,c and Xt are given by (|23|) . (|20p and (fT3| respectively. 

Now let u'^{t, X, y),u{t, x) be the solutions to ([T2l) and (|19p respectively. Tak- 
ing into account relations ([25)1 . (1^^ . the weak convergence of X^ to Xt as e — > 
(which is implied by Lemma 3.3) and Proposition 3.2 we get the statement of 
the Theorem. 

n 

We conclude this section with the a-priori bounds for the Holder norm of the 
solution and for the sup-norm of the solution, the first and the second derivatives 
of the solution of ([T]). These bounds will be uniform in e. The method follows 
closely [S]. 

Let us first introduce some notation. 

We write U^ = [0,T) x D% t/^ = [0,T) x D% dU^ = [0,T) x dD' and 
Vf = (0, T)xD\ where U' = Z?' U dD\ 

For 0<a<l,T>0 and for any function g we write: 

||5llc/|, = sup \g{t,z)\ 
\\H^9\\u, - sup l^(M)-5(t',^')l 



(M),(t'"')ec/» ii-tr/' + k-^r 
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\\g\\u^,a = Mu^ + WH'^gWu^ 

\\g\\D',T,l+a = \\9\\u:j,.a + \\gt\\ui. + \\Dg\\(o,T)xD^ 
ll5llD%T,l+a = \\g\\D',T,l+a 

Moreover for notational convenience we will write z = {x,y). 

Lemma 3.5. Under our assumptions there exists a constant Ci, indepen- 
dent of e > 0, such that 

< u' < Ci in Ut. 

Proof. Lemma 3.5 can be proven using equation (J12p . Here we give an 
analytic proof of the claim. For any fixed 6 > we define the function 

w" = {u" - b)+ = max{u' - b, 0}. 

It is easy to show that 

wl < -Aw' on (0,T) X D\ 

in the weak sense. Let us choose now b = max{A^, ||/||}, where N is such that 
if u > TV then c(-, •, u) < -M for some M > 0. Then 

w'iO,x,y)=0. 

Let us now assume that w'^ attains a maximum positive value on the bound- 
ary dV^ at the point {to,Xo,yo)- Since u'^ is continuous up to the boundary, 
there exists a connected set A such that {to,Xo,yo) S A, A C dV,^ and w"^ > 
on A, i.e. w*^ > 6 on A. Since {to,Xo,yo) is a maximum for w'^ and 7'^ is the 
inward normal derivative we get that 4^ < at {to,Xo,yo)- But on A we have 
that ^TT = f^ = —£c{x,y,u'^)u'^. Taking into account the particular choice 
of b and that c(-,-,u) < —M for u > N, we get that —ec{x,y,u'^)u'^ > at 
(to,Xo,yo)- Thus we have a contradiction and so maximum principle implies 
that 

w^ = => u" <b in [/^. 

Lastly maximum principle again implies that m*^ > 0. 

n 

Let us consider the following linear parabolic pde: 

vl = ]-Av% m{0,T)xD' (27) 

«'(0,x,2/) = fix), on{0}xl?^ 

-^ = ^ec{x,y)v\ on{0,T)xdD', 

10 
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where /, c are bounded smooth functions. Under the standard hypotheses prob- 
lem ((27)) has a unique classical solution (Theorem 5.3.2 in [9]). 

Lemma 3.6. There is a constant C, independent of e, and an open set 
/ C (0, 1) such that for any a £ I: 



M\D',T,l+a + \\D'vlv^<C. (28) 

Proof. We will give just a sketch of the proof, since the analysis follows [14], 
[IS] . [TO] and [5]. The calculations are lengthy but standard. 

We solve the second initial-boundary value problem (1271) by reducing it to 
an integral equation, i.e. we write: 

v'{t,z)^ f f T'{t,z,T,0<f>'{T,OddDldT+ [ r'{t,z,0,Ofmi, (29) 

Jo JdD' JD' 

where r<=(i,z,r,0 = (2y^)-"-"(i - t)-"^ exp[--^2i^^^^2_ii)] is the funda- 
mental solution to the heat equation and ^(i, z) is the solution to a Voltera type 
integral equation: 



JdD- dl ^ 

Jd- oY Jd' 



Let us now define 



F\t,z) = / ^EllMi^f(^^)d^ + ,c{z) [ r(i,z,o,o/(C)rf^ 

Jd- (^1 J D' 

M,{t,z,T,0 = ^£!Mlli)+,c(z)r(t,z,r,e) 

M^+i(t,z,T,<e) = f f Mi{t,z,t',z')M^{t',z',T,^)ddDl,dt' 

Jo JdD- 

It can be shown (see f5]) that there is a Holder continuous (in space variables) 
and bounded (with bound and Holder coefficient independent of e) solution cj)'^ 
for (|30|. expressed in the form: 

OO n^ n 

(j)''{t,x)=2F'(t,z) + 2y2 / M^{t,z,T,C)F'{T,^)ddDldT (31) 

^^iJo JdD' 

Using the boundedness and the Holder continuity of (PT|) and (|29l) , one can show 
(see [M], [15], [16] and [9]) that there is a constant C, independent of e, such 
that 



lD^T,l+o 



l^'^^'lk^ < C. 



11 
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n 



Now we are ready to prove the result for the a-priori bounds: 



Proposition 3.7. There is a constant C, mdependent of e, and an open set 
/ C (0, 1) such that for any b > a E I {a is the constant from Lemma 3.6.): 



\\uqn'-.Ts+b + \\D'u^H<C. (32) 

where u' is a classical solution to ([T]). 

Proof. We will use Schauder's fixed point Theorem. Let us first define for 
convenience || • ||2_^^ = || • ||£„_r_i+„ + \\D^ ■ \\vf- 

Let C^+a ^g ^^^ Banach space of all functions u'^{t, z) that are continuous in 
Urp with norm Hu'^Hj i „• 



For any C > 0, let 6^+" be the set {u' : u" G 6^+", \\u%^^<C}. 
For every u' S 6^'*''' define w' = Tu^ to be the solution to the following 
problem: 



wl = |aw% in(0,r)x£)'^ (33) 

w\Q,z) ^ fix), on{0}xi:>" 

-— = -ec(z,u<^)w% on (0,T) X aL*', 



^7' 
Then similarly as in Lemma 3.6, one can write: 



w'it,z)= [ f V'{t,z,T,iW{T,OddDldT+ [ T%t,z,0,Of{Odi, (34) 

Jo JdD' JD' 

where 4)'^{t,z) satisfies: 

(p\t,z) = 2f f [^Elihlfl^ + ec{z,u^)r^{t,z,T,OW{T,OddDldT 
Jo JdD- dl * 

+ 2[/ S|£lM)/(^)d^ + ec(z,u^) / r{t,z,0,0f{0dm5) 

JD- CI JD' 

We shall prove that T has a fixed point. 

Since u*^ and c are bounded functions, one can show, in the same way as in 
the proof of Lemma 3.6, that the function 4f{t, z) that satisfies ([35]) is bounded 
and Holder continuous (in space variables) with bound and Holder constant 
independent of e. 



12 
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Using this result and representation p4p one can conclude (Lemma 3.6) that 
there is a constant C such that 



\\w%+a < C. 

So T maps Cp into itself for an appropriately chosen constant C. 



Now let {u^} be a sequence of functions that belong to Cq'^ and 'w^,4>%^ 
be defined by pi]) and p5|) when u^ = u^. Assume that ||u^ -- "^Il2+a ^ as 
n — >■ oo. We need to show that ||w^ — '"^'^112+a — >■ as n — ?► cx). 

The continuity of the function c{z, u) in w- variables imply that || '/'^— </>' || c/" — >■ 
as n — 7* oo. This and ([M|) give us ||z«^ — ^"^112+0 ~^ 0. 

Therefore T is a continuous map. 

Next we need to show that T maps C^"^" into a compact subset of Cp"*"". 
This is an easy consequence of Theorem 7.1.1 of [9], which states that for < 
a < b < 1, the bounded subsets of C^+'' are pre-compact subsets of C^+°. 

Lastly G"^ is a closed convex set of the Banach space 6^+''. 

Therefore by Schauder's Fixed Point Theorem we get that T has a fixed 
point, i.e. there exists a u'^ such that u"^ — Tu'^ and actually 

n 



4 Some Results On Wave Front Propagation 

In this section we will see some applications of Theorem 3.4 to the question of 
wave front propagation in narrow domains. As we mentioned in the introduc- 
tion, corresponding results on the standard reaction diffusion equation ([2]) (see 
chapter 6 and 7 in [4], [10] and [13]) allow to describe the asymptotic wavefront 
motion for H]). 

We will focus on two different cases. In subsection 4.1 we consider the 
case where the functions c(-,u), V{-), S{-) and /(•) change slowly in x, i.e. 
c(.,u) = c{5x,u), V{-) = V{5x), S{-) = S{5x) and /(•) = f{5x) for < (5 < 1. 
We first assume that the nonlinear boundary term in ([1]), c{x,y,u), is of K- 
P-P type for y = 0, i.e. c{x,0,u) is positive for u < 1, negative for u > 1 
and c{x) = c(a;,0,0) = maxo<u<i c(a::, 0, u). We will see how the motion of 
the wavefront depends on the behavior of the cross-sections D^ of the domain 
D. In particular, using the results of 4| (chapter 6) we will see that in the 
case of the nonlinear term of K-P-P type and for a; G M the wavefront can 
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have jumps. Actually, the jumps of the wavefront appear at positions where 
the tube becomes thinner. The results are given in Theorem 4.1.2, Theorem 
4.1.5 and Theorem 4.1.7. Then we briefly discuss the bistable case, i.e. when 
c{x,0,u) > for li G (/i, 1) and c{x,0,u) < for u e (0,/^) U (l,oo), where 
< /i < 1. In this case we consider a specific example and we will see how 
the asymptotic speed of the wavefront depends on the surface area to volume 
ratio T77> ■ In subsection 4.2, we return to the K-P-P case, but now we consider 

V(x) ' ' 

front propagation when x S M and the boundary dD^ of D^ is determined by 
stationary random processes on M on some probability space (0,5", P). The 
conclusion is in Theorem 4.2.7. 

We will denote by c{x,u) := 2v&:c{x,Q,u{t,x)) the nonlinear term in ([2]). 
Obviously the type of c(a;, u) (K-P-P or bistable) is determined by c{x, 0, u) and 
vice- versa. 



4.1 Wave Fronts in Slowly Changing Media 

Let us assume that the functions c{-,u), V{-), S{-) and /(•) change slowly in 
x, i.e. c{-,u) = c{dx,u), V{-) = V{Sx), S{-) = S{6x) and /(•) = f{Sx) for 
0<(5< 1. 

We start with the case where the nonlinear term c{x, u) of ([2|) is of K- 
P-P type. We additionally assume that the closure of the support of /, Fq, 
coincides with the closure of its interior. Lastly we take for brevity a; e M^ and 
c(x) = c(x, 0) — ^Yn^c{x,0,0) (recall that c(x,0, 0) = supQ^^^ic{x,0,u)) to 
be an increasing function. 

Let (p : [0, T] — >■ R^ and introduce the functional 

D fj,\ ) L [c{4>s) - h\'Ps\'^]ds, (/> is absolutely continuous . . 

I +00, tor the rest ot Cq.t- 

Put 

W{t, x) = sup{i?o,t(0) : e eo,i(Ki), 00 = x, 0* e FJ. (37) 

We say that condition (N) is satisfied if for any t > and {t,x) E {{t,x) : 
W{t,x) =0} : 



W{t, x) = sup{i?o,t(0) : 00 = X, <i)t e Fo, (i - s, 0,) e {(i, x) : W{t, x) < 0}}. 

As it is mentioned in chapter 10 of [8[, condition (N) is fulfilled for the 
smooth and increasing function c{x). Moreover as we shall see in Theorem 
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4.1.2, W{t, x) determines the motion of the wave front for u'^ for small enough 
e>0. 

Let us consider u{t,x), the solution to equation ^, for n = 1. If we set 
u^{t, x) = u{t/6, x/6), then u^ is the solution to the following parabolic problem: 

uf = ttuL + o T77-r"!^ + tc{x, u^{t, x))u\ in (0, oo) x R^ 
2 2 V{x) 5 

u\0,x) = f{x) > 0, on {0} X R\ (38) 

Under the assumptions above, the following theorem, which is a reformula- 
tion of Theorem 6.2.1 of [4], states that W{t,x) determines the motion of the 
wave front for u {t,x) under condition (N): 

Theorem 4.1.1. Let u^{t,x) be the solution to ([38|) . Then under condition 
(N) we have: 

sio ^ ' ' \0, W{t,x) <0. ^ ' 

Let us consider now equation ([T]) for n = 1, c(-,m) = c{Sx,u), /(•) = f{Sx) 
in a slowly changing in x narrow domain D*^'^ , so that V{-) — V{Sx), S{-) = 
S{Sx). Let us define u'^'^{t, x, y) = u'^{t/6, x/S, y). Under the assumptions above, 
Theorems 3.4 and 4.1.1 imply that W{t, x) will determine the motion of the wave 
front in this case too, as follows: 

Theorem 4.1.2. The following statement holds: 

limlim.^-^(M,y) = r' ^(^'-)>0 (40) 

SiO eiO ^ ' '^^ \0, W{t,x) <0. ^ ' 

So the equation W{t, a;) = defines the position of the interface (wavefront) 
between areas where u^''^ (for e > and (5 > small enough) is close to and to 
1. Actually, as we shall see below the wavefront may have jumps. It is known 
(see chapter 6 in [J]), that because of the dependance of c{x) on x, the wave 
front of u may have jumps and new sources may be "igniting" ahead of the 
front. We will give sufficient conditions that guarantee such jumps for a class 
of smooth and increasing functions c{x). Hence Theorem 4.1.2 implies that one 
can predict appearances of new sources and jumps of the wave front of w'^' for 
e > and (5 > small enough. Reaction-Diffusion Eqautions (RDE's) with 
Nonlinear Boundary Conditions In Narrow Domains Let t* = t* (x, c(-)) be such 
that W{t* ,x) = 0. Such a t*{x,c{-)) is defined in a unique way. 
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Figure 1: A wavefront that jumps from xq to X2 at time to. 

We have the foUowing proposition (see chapter 6 in [4] for more details): 

Proposition 4.1.3. Let t*{x) be as in Figure 1 and Fo — {x e R^,x < 0}. 
Then the wavefront jumps from Xo to X2 at time to (see Figure 1), i.e.: 

(i). li t < to then hm^io linie^o u'' (t, x,y) ^ 1 for a connected set: 
Ft = {x £R^ : W{t, x) > and a; < xq}. 

(ii). If <o < ^ < ^1 then the set where hm^^o linic^o w*^' (i, x, y) — 1 consists of 
two connected components: 
Ft = {x e M^ : W{t, x) >0 and x < xi} U {x e M^ : W{t, x) > and x > 

Xl}. 

The set {x € M^ : W{t,x) > and x < xi} is at a positive distance from 
the set {x e K^ : W{t, x) > and x > xi} for to < t < ti. 

(iii). If i > ti then Hm^^o lime^o u'^ {t, x,y) = 1 for a connected set: 
Ft = {xeKi :W{t,x) >0}. 

Based now on comparison resuhs (Lemma 4.1.4) we wih give sufficient con- 
ditions that guarantee jumps of the wavefront. In particular we will prove 
(Theorem 4.1.5) that if c(x) is a rapidly increasing smooth function, then 
t* = t*(x, c(-)) such that W{t*,x) = is as in Figure 1. 

The functional Ro.xi'P) defined in ([M)) and the function W(t,x) defined 
in (I37|) depend also on c. Hence we will write sometimes i?o,T(07 c(-)) and 
W{t,x,c{-)) in order to emphasize this dependence. 

We have the following comparison result: 

Lemma 4.1.4. 
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(i). Let A be a positive number. Then t*{x,Ac{-)) = -j=t* {x , c{-)) . 

(ii). Let a be a positive number and define Ca{x) — c{ax). Then t*(x,Ca{-)) = 
H*{ax,c{-)). 

(iii). Let ci, C2 be two functions such that ci{x) < C2{x) for every a: £ M^. Then 

t*{x,Ci{-))>t*{x,C2{-)). 

Proof. Let us write t^ = t*{x,Ac{-)) and let (j)^ be the extremal so that 
W{t\, X, Ac{-)) = Ro,t*^ (0'^, Ac{-)) = 0. Such an extremal satisfies the following 
Euler-Lagrance equation: 

4>^(s) - -Ac'icl^^is)) 

(f>'^{Q) = X (41) 

^^it*A) - 0. 



Let us define now the function 0(s) — (p"^ {s / ^/A) . We claim that the function 
(j}{s) is the extremal so that W{^/At*y^,x,c{-)) — i?o V3t* i4':C{-)) = 0. Indeed 
it is easy to see that the definition of (p and the fact that i?o.t* {4>^, Ac{-)) = 
imply that i?„ ^/At' (^^ ^i')) = 0- Moreover (/) satisfy an Euler-Lagrange equation 
of the form (HT|) with Ac{x) and t\ replaced by c{x) and VAt\ respectively. This 
proves the claim, which implies part (i) of the lemma. 

Part (ii) of the lemma can be proven in a similar way. We define t* = 
t*{x,Ca{-)) and let (/)° to be the extremal so that W{t'^,x,Ca{-j) — i?o,t* (0", Ca(-)) = 
0. Then similarly as it is done in part (i), one should consider the function (p{s) 
that is defined by 0(s) = a4>°-{s/a). 

We prove now part (iii) of the lemma. Let us define t\ = i*(x,ci(-)) and 
i* = t*{x,C2{-))- Moreover let 0^ be the extremal so that W{t\,x,ci{-)) = 
Ro,ti{<P^,ci{-)) — 0. Since ci(x) < C2{x) we have 

- Ro,ti (0\ ci(-)) < Ro.ti (0\ C2(-))- (42) 

Furthermore, it is easy to see that W{t, x) is an increasing function of t. 

Let us assume now that tl < ^2- This assumption and the fact that W{t2,x, C2(-)) 
imply that W{tl,x, C2(-)) < 0. By recalling the definition of function W, one 
easily concludes that: 

i?o,n(0\c2(-))<O. (43) 

However inequality (|43|) contradicts (|42|) . Therefore t*{x,ci{-)) > t*{x,C2{-)). 

D 
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In section 6.2 of fT, it is proven that if c{x), instead of the smooth function 
2 y7^c(a;, 0, 0), is a piecewise constant function, denoted by d{x), such that 

Jf \ ( di, X < X2 

d{x) ^\ (44) 

[d2, X > X2. 

with(i2 > 2di > 0, then the function i* = t*{x,d{-)) such that W{t*,x,d{-)) = 
is not monotone, as in Figure 1. More specifically the curves connecting the 
point (0,0) with {xi,ti) and (xi,ii) with (0:2, io) are line segments and for 
X > X2, t* — t*{x, d{-)) is the solution to 

sup{d2(r -t) + d,t %,^ - il = 0. 

Moreover in this case 



2{t* - t) 2t 



V2{d2-di) 

to = X2 1 (45) 

"2 

ti = --—={x2 + ^2d'ito) (46) 

Z\' zdi 

We will write to = io('^) and ti == ti{d) to emphasize the dependence of to and 
ii on the function d(x). 

With the help of the result above and Lemma 4.1.4 we will give sufficient 
conditions that guarantee jumps of the wavefront of u {t,x) (and by Theorem 
4.1.2 of u^(t, X, y) for e > and S > small enough) for a class of smooth and 
increasing functions c. 

Let us define the set 



A = {{di,d2) eR\xR\:d2> 2di and ^2 > 2^/diid2 ~ di)}. (47) 

It is easy to see that A is a non-empty set. 

Theorem 4.1.5. Let d{x) be the step function defined in (^1)) such that 
{di,d2) G A. Consider real numbers A and a such that 

(i). a,A> 1. 

(ii). aVA < i[l + ^ =]. 

Then for any smoothly increasing function c{x) such that 

d{x) < c{x) < Ad{ax) (48) 

the wavefront corresponding to c has jumps. In particular the excitation reaches 
the region {a; > —+6} before it reaches the point — , where 6 is a. small enough 
positive number and — is as in Figures 2 and 3. 
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Proof. Let us define d{x) = Ad{ax). Since a,A > 1, the function d(x) is 
shifted vertically upwards and horizontally to the left. So we get that d{x) < 
d{x) (see Figure 2). 

Parts (i) and (n) of Lemma 4.1.4 imply that t*{x,d{-)) = —j=t* {ax , d{-)) . 
This and part (iii) of Lemma 4.1.4 give that if c satisfies (|3S|), then t*{x^c{-)) 
will satisfy (see Figure 3): 



a\/A 



t*{ax,d{-)) < t*{x,c{-)) < t*{x,d{-)). 



(49) 



FLjnc^tion c^r) 



Position of the w3vefront 







Figure 2: d{x) < c{x) < Ad{ax) Figure 3: t*{x,d) < t*{x,c) < t*{x,d) 

We know that t*{x,d{-)) and t*{x,d{-)) are not monotone (recall that d and 
d are piecewise constant functions). We will show that t*{x,c{-)) is also not 
monotone (i.e it is as in Figure 1). Let us assume that 



h{d) > to{d), 



(50) 



where tQ{d) is as in psj) and ti{d) is defined similarly to ti{d) in (H51) with 
di , ^2 , X2 replaced by Adi , Ad2 , ^ respectively. In particular ([50]) holds if con- 
dition (ii) above holds, i.e. if a\/A < ^[1 



2^ 2y^di(d2-di)' 



Moreover the con- 



dition ^2 > '2\/di{d2 — di) in the definition of A in (|47p implies that i[l 



d2 



> 1, which has to be true since a, A > 1. 



2y/di{d2-di)' 

Inequality ([50]) can be equivalently written as i*( — ,d(-)) > t{x2,d{-)). By 
this and (HHl) we immediately get that 



t*{x2,c{-))<t*{^,c{-)) 

a 



(51) 
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which, since — < xi < X2, implies that t*{x,c{-)) is as in Figure 1 and so new 
sources are igniting ahead of the wave front. 

In Figures 2 and 3 we see an illustration of the construction. 

n 

Example. An example of a function c{x) that satisfies the requirements of 
Theorem 4.1.5 is 

cix) — TT — n , (52) 

where (di,(i2) G A, a, A satisfy assumptions (i) and (ii) of Theorem 4.1.5, 
k G { — ,X2) and the constants ^ and A are chosen so that c( — ) < Adi and 
c{x2) > d2- 

In particular now if c(x) = ^ yfl , i.e. c(a;, 0,0) — 1, is an increasing smooth 
function that satisfies the requirements of Theorem 4.1.5, then the jump of the 
wavefront of u'^'^{t,x,y), for e > and 6 > small enough, occurs when yfi 
increases rapidly. This implies, at least when the tube D^ retains its shape as x 
increases, that the jumps of the wave front occur at places where the tube D^ 
becomes thinner, i.e. when V{x) decreases significantly. 

Remark 4.1.6. Similar results hold for layers as well, i.e. for x E K" with 
n > 1. 

Using the results in f7] one can consider the limiting behavior as 5, e J, of 
u^'^{t, X, y) when condition (N) is not fulfilled. We will briefly discuss the result 
for the general case x € M". 

Instead now of function W{t,x) defined by (|37l) . we consider the function 

W*{t,x) = sup{ min i?o,s('/') : S Co,t(]R") is absolutely continuous, 

0<s<t 

00 = a;, (t>teFo}. (53) 

One can prove that W* {t, x) is Lipschitz continuous and that W* {t, x) < 
mm{0, W{t,x)}. 

Then Theorem 2.1 in [7] and Theorem 3.4 imply that W*{t,x) determines 
the motion of the wave front as follows: 

Theorem 4.1.7. The following statements hold: 
(i). For any compact subset 0i of the interior of {{t, x) : t > 0, W*{t, x) — 0}, 
limlimM*^' {t,x,y) = 1 uniformly in {t,x) £ 9i. 
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(ii). For any compact subset 62 of {(t,x) : W*(t,x) < 0}, 

limlimu'^' {t,x,y) = uniformly in {t,x) G 82- 
5^0 €^0 



We conclude subsection 4.1 with the case that the nonlinear term c{x,u) 
of ^ is of bistable type, i.e. c{x,u) > for u e (m^I): c{x,u) < for w € 
(0,/i) U (1,00), where < /i < 1. This problem was considered in [10] and it 
was also presented in section 6.4 of jl]. 

Here we restrict the analysis to a concrete example that allows to give an 
exact formula for the asymptotic speed of the wavefront of u'^'^ for e > and 
(5 > small enough. As we will see the asymptotic speed of the wa vefro nt is 
proportional to the square root of the surface area to volume ratio \ yf\- 

To be specific let x € M", c{x, 0, u) = (u — /z)(l — m), < ^ < ^ and assume 
that the function u^{t,x) (compare with ([38])) is the solution to 

-^ ~ ^ -div{V{x)V^u^) + ]-^^iu^ -fi)il-u^)u\ in(0,oo)x 



' 2V{x) ^ ^ ' - ' 6 2Vix) 

u\0,x) = f{x), on{0}xR". (54) 

Consider a point x G M" to be excited at time t, if u^{t,x) (the solution to 
([M]) ) is close to 1 and non-excited if u*(i, x) is close to 0. Then the Corollary of 
Theorem 4.1 of [10] gives us that for small S > Q the region {x e M" : f(x) > fi} 
becomes excited and the region {x £ K" : f{x) < /i} becomes non-excited after 
a short starting phase. Now let u'^'^{t,x,y) = u'^{t/6,x/S,y), where u'^{t,x,y) 
is the solution to ([T]). Theorem 3.4 implies that the same conclusions hold for 
u^'^{t, X, y) for e > and ^ > small enough. 

To compute the asymptotic propagation speed of excitation at a; G M", let 
us consider the equation for the wave profile: 

\vl^{0 + a{x)v'^{i) + ^||y(«(0 - M)(l - v{0)v{i) = 0, C e ^(55) 
lim i;(^) — 1, lim v{S,) = 0. 

As it can be verified by direct substitution, equation (|55|) is solvable if a{x) is 
given by the formula 

Moreover, in our case, (15^ is also the asymptotic propagation speed of excitation 
at a; € M" and it is independent of direction. 
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Lastly, it is known that as the size oi D^ increases (without changing shape), 
the surface area to volume ratio y0: decreases. In the case x € R, this fact, 
equation ([5S|) and Theorem 3.4 imply that the wavefront of u"^'^ (for e > and 
(5 > small enough) slows down when the tube becomes thicker. A similar 
result holds for layers. 



4.2 K-P-P Fronts in Random Media 

In this subsection we consider wave front propagation for the solution of ([T]) for 
small e > 0, when x G M, the boundary dD^ of D^ is determined by stationary 
and ergodic random processes on M and the nonlinear boundary term in ([T]) (for 
y = 0, i.e. c{x,Q,u)) is of K-P-P type. As we did in subsections 4.1, we will 
first consider (Theorem 4.2.6) wavefront propagation for the solution of ([2|) and 
then with the aid of Theorem 3.4 we will consider (Theorem 4.2.7) wavefront 
propagation for the solution of ([T|) for small enough e > 0. As we will see the 
cross sections Dx of D affect the speed of the wavefront through the surface to 
volume ratio Ytx)- 

In sections 7.4 — 7.6 of [4] wave front propagation for equations like Q is 
considered in the case where there is no drift term and the randomness comes 
only from the nonlinear part of the equation. Moreover in |13j the authors 
considered the case of reaction-diffusion equations of type ([2|) with a random 
drift and homogeneous in x nonlinear term. In the case considered here, both 
the drift and and the nonlinear term are random. In [3], pp. 524-525, the author 
remarks that one could use the procedure developed in sections 7.4 — 7.6 of [4] 
to study wavcfronts in one-dimensional uniformly bounded random drift with 
random nonlinear term. We will see that one can prove Theorem 4.2.6, which 
is analogous to Theorem 7.6.1 in [i], by following the proof of Theorem 7.6.1 
in [4] . We make use of the results in [13] and of the fact that the operator of 
the equation ([2]) is self adjoint with respect to an appropriate inner product 
(it has the form 2^/Tx)l^(^i•'^)l^))■ Actually the latter simplifies the analysis 
significantly. Instead of repeating the proof of [4 here, we will only outline the 
differences. 

Let us first list our assumptions. Consider a probability space (17,5^, P). 
We assume that the random field V{x,(li) (namely the volume) is three times 
continuously differentiable, i.e. V £ C'^(R), with P probability one. Suppose 
that Q{x) = {■£:(\ogV{x)), y^) is a random vector function on {Cl,^,P) and 
that it is measurable, stationary in x and translation in x generates an ergodic 
transformation of the space fl. Moreover the function ■^{\ogV{x)) is assumed 
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bounded, with zero mean (i.e. E[-^{\ogV{x))] = 0). We additionally assume 
(for the purposes of Lemma 4.2.1 and 4.2.3) that there is a set of nonzero P 
probability on which 

lim I [V{x,io)]-'^dx = +oo. (57) 



If condition (j57l) holds on a set of nonzero measure then, by the ergodicity 
assumption, it must hold with P probability one. 

As far as the non-linear term c(x,u,lu)u = yf'^l c{x,0,u)u is concerned, 
in addition to the stationarity and ergodicity assumptions, we also make the 
following assumptions. For all x S M, c is of K.P.P type, i.e. c{x,0,u) is 
positive for u < 1, negative for m > 1, continuous in u for u > and c{x) = 
c(x, 0,0) = supq^„ c(x, 0, u). Moreover with P probability one, the function 
c(a;, u, Cj)u satisfies a Lipschitz condition of the form 

\c{x,Ui,Uj)ui -c(x,U2tUj)u2\ < ^^, ' . C,{x)\ui - M2I, for X,Ui,U2 G K, 

V{x,uj) 
such that for alH > and a; e M, 

E-^mI |7^C(^.)} < oo, P-a.s., 
Jo y(^s) 

where {Xt,Px) is a diffusion process with random generator L ~ ^tt + 
l£ilogVix,u.))£. 

The initial function f{x) is assumed to be nonnegative, bounded from above 
and non-random. 

Let now ^{z) be the function defined by the equality 

fi{z) = ^[lni;ixro<oo exp{ / ' [ciXs) + z]ds}], z e M, (58) 

Jo 

where c{x) = ^TFKr)C{x, 0, 0) and tq is the first hitting time of the process Xt to 
the point 0. For tq one has the following lemma: 

Lemma 4.2.1. Condition dST]) and E[^{\ogV{x,uj))] = imply that 
Pi(ro < oo) = 1. 



Proof. It follows directly from the proof of Lemma 4.4 of [13 if one notes 

1 _d_/ 

2 dx^ 

n 



that the drift term is 5-^(logy(a;)) 



Lemma 4.2.2. Under the assumptions imposed above, function /i(z) has 
the following properties: 
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(i). For all z € R, ^(z) = limj^-oo 7 ln£'fXTo<oo exp{/J"°[c(Xs) + z]ds}. 

(ii). Function fi{z) is convex, lower semicontinuous and monotonically non- 
decreasing in z. Moreover fi{z) is continuously differentiable and the 
derivative n'{z) is positive and monotonically increasing for z < g^, where 
(jf^ is a non-positive number (which actually is the discontinuity point of 
/i(z), as property (iii) below shows). 

(iii). ii{z) < for z < g^ and ii{z) = 00 for z > g^ where 5^1 < 0. 

Proof. Property (i) can be proven as Proposition 2.1 of [T3]. Property (ii) 
follows similarly as Theorem 7.5.1(ii) of 4 . Property (iii) follows analogously 
to Theorem 7. 5.1 (iii) of [1]. Here one uses the fact that the operator of ^ has 
the form 2i7(^^C'^(^)£))' i-^- it is self adjoint. 

n 

We also observe that /^(z) > ^o{z) where Hoiz) — E[\n Ei{xTo<ocie^'^°)]- As 
it has been proven in Lemma 2.2 of |13| . function ijLo{z) has properties (i)-(iii) 
of Lemma 4.2.2 as well (for c{x) = 0). In addition the following lemma holds, 
which is a restatement of Proposition 4.1 of J13j . 

Lemma 4.2.3. Condition dST]) and E[-^{\ogV {x , uj))] = imply that the 
discontinuity point of /io(z) is 5^„ = 0. 

We will assume that —00 < g^ <0 (by Lemma 4.2.2(iii) or Lemma 4.2.3 we 
already know that 5^ < 0) and we define I{y) = sup2<^ [yz — //(z)] for ?/ G M. 

Lemma 4.2.2 and the fact that ^i{z) > ^o{z) imply that the arguments in 
the beginning of section 7.6 of 4J carry out here as well. Therefore we conclude 
that there is a unique v* > such that /(-^^r) = and i/* = inf2<g^_ 'Tz)' 

Remark 4.2.4. We would like to emphasize that the existence and unique- 
ness of a positive ly* follows mainly from properties (i)-(iii) of fj,{z) (Lemma 
4.2.2). In particular property (iii) holds because the operator of ([2]) is self ad- 
joint. 

Similarly as Theorem 7.6.1 in [1] was proven, one can prove Theorem 4.2.6 
below. 

Note that by following the proof of Theorem 7.6.1 in [4], one needs to esti- 
mate certain probabilities for tq and Xt . For this purpose we have the following 
lemma: 

Lemma 4.2.5. Let 5 be a positive number and Us{0) = {x : \x\ ^ 6}. Then 
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(i). inf^e[/,(o) Px{to s^ 1} > 0, P-a.s. 
(ii). infa;e[/,(o),se(oa] Px{Xs e Us{0)} > 0, P-a.s. 
(iii). For a > and rj > S > we have 

inf PAr-n-a > 1,^1 e Us{0)} > 0, P-a.s. 

xSiUs( — a) 

Proof. The proof of all statements follows from the corresponding state- 
ments for Wf in place of Xt (see for example section 7.5 of ^) and by the 
Girsanov's theorem on the absolute continuous change of measures in the space 
of trajectories. 

n 

Therefore we have the following Theorem: 

Theorem 4.2.6. Let x e M and u{t,x) satisfy equation ([5]). Under our 
assumptions we have: 

(i) . For all i^ > ly* , 

lim sup u{t, x) = 0, P — a.s. 

*->-o°£c>jyt 

(ii). Let us define Ch{x) = ^yf^ info<M</i c{x, 0, u) and assume that there is a 
constant k> such that for any < h < 1 and a; G K, 

K < Ch{x), P — a.s. 

Then for all t^e (0,j^*), 

lim sup u{t,x) — l, P ~ a.s. 

Finally Theorem 3.4 and Theorem 4.2.6 imply: 

Theorem 4.2.7. Let {x,y) e M x M™ and u''{t,x,y) satisfy equation ([T|). 
Under our assumptions we have: 

(i). For all i^ > ly* , 

lim sup lim u'^(t, a:, y) = 0, P — a.s. 
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(ii). Let us define Ch{x) = ^jrt§) info<u<h c{x, 0, u) and assume that there is a 
constant k > such that for any < h < 1 and a; G M, 

K < Ch{x), P — a.s. 

Then for all i^ & {0,v*), 

lim sup lim u^(t, a;, y) = 1, P — a.s. 

t-i-oo Q<2,<^( e-yO 

Remark 4.2.8. Theorem 4.2.6 was proven in ([4]) with the assumption in 
part (ii) replaced by the assumption that for any < h < 1 and i/ e M, 

1 /•* 

limsup -lni?i,t exp{— / Ch{Xs)ds} < 0, P — a.s, (59) 

t— >00 t Jq 

which is however difficult to verify. Obviously the assumption made in part (ii) 
of Theorems 4.2.6 and 4.2.7 implies §9i). 
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